Tracing light propagation to the intrinsic accuracy of space-time geometry 
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Advancement in astronomical observations and technical instrumentation requires coding light 
propagation at high level of precision; this could open a new detection window of many subtle rela- 
tivistic effects suffered by light while it is propagating and entangled in the physical measurements. 
Light propagation and its subsequent detection should indeed be conceived in a fully relativistic 
context, in order to interpret the results of the observations in accordance with the geometrical 
environment affecting light propagation itself, as an unicum surrounding universe. One of the most 
intriguing aspects is the boost towards the development of highly accurate models able to recon- 
struct the light path consistently with General Relativity and the precepts of measurements. This 
paper deals with the complexity of such a topic by showing how the geometrical framework of models 
like RAMOD, initially developed for astrometric observations, constitutes an appropriate physical 
environment for back tracing a light ray conforming to the intrinsic accuracy of space-time. This 
article discusses the reasons why RAMOD stands out among the existent approaches applied to the 
light propagation problem and provides a proof of its capability in recasting recent literature cases. 
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I. INTRODUCTION 

The treatment of light propagation in time-dependent 
gravitational fields, as the surrounding universe, is ex- 
tremely important for astrophysics, encompassing is- 
sues from fundamental astronomy to cosmology. Attain- 
ing very accurate measurements will allow to observe a 
new range of subtle physical effects (see, for example, 
Will [H, Turyshev 0, Crosta and Mignard 0, Kopeikin 
and Gwinn ML Kopeikin et al. p, |6|, Blanchet et al. 
, Kopeik in p| , Kopeikin and Makarov Q , Damour and 
Nordtvedt pTol ]T Damour et al. (Tlj)- 

Undoubtedly, light tracing should be conceived in a 
general relativistic framework. Today General Relativity 
(GR) is the theory in which geometry and physics are 
joined together in order to explain how gravity works. 
The trajectory of a photon is traced by solving the rela- 
tivistic equations of the null geodesic in a curved space- 
time and, at the same time, the detection process usu- 
ally takes place in a geometrical environment generated 
by a n-body distribution as could be that of our Solar 
System, i.e., by a localized gravitationally-bound distri- 
bution of masses moving with their velocities and accel- 
erations. Once the physical set-up is defined, the null 
geodesic should be made explicit according to the corre- 
sponding geometry. 

Nevertheless the exact solution of the null geodesic is 
not generally known. Hence, one has to resort to ap- 
proximation methods. Nowadays several conceptual ap- 
proaches model the light propagation problem in the rel- 
ativistic context. Among them, the post-Newtonian (pN) 
and the post-Minkowskian (pM) approximations are 
those mainly used (Kopeikin and Schafer [12], Kopeikin 
and Mashhoon fljj . Klioner [ljj], Teyssandier and Lc 
Poncin-Lafitte [15j : and references therein). 

The peculiar between the two approximations consists 
in the fact that pN is based on the assumption of an 



everywhere-weak gravitational field and of slow motion, 
while the pM only on the weakness of the gravitational 
field. This can be formalized [Hj] by saying that the 
expansion of the metric tensor can be characterized in 
powers of a small parameter, 1/c in the pN case [ItJ and 
G in the pM approximation (which contains also terms 
in 1/c). Basically the approximation methods rely upon 
the idea of being pN as close as possible to the New- 
tonian theory (absolute time, absolute space, auxiliary 
Euclidean metric and instantaneous potential) and pM 
to the Minkowskian interpretation of special relativity 
(absolute space-time, auxiliary Minkowskian metric and 
retarded potential) fig . 

In the framework of the pN approximation of GR, the 
light path is solved using a matching technique which 
links the perturbed internal solution inside the near- 
zone of the n-body system with the external one. Ap- 
plications of this method are adopted, for example, in 
HH EM EH j where the boundary conditions for the null 
geodesic are fixed by the coordinate position of the ob- 
server and imposing the value of c to the modulus of 
the light ray vector at past null infinity. However, the 
expansion of the metric tensor may not be analytic in 
higher pN approximations (l6| . The second pN approxi- 
mation in propagation of light was investigated by Brum- 
berg in [20(. Advancement on second order pN equa- 
tions of light propagation are made in Xu et al. [2l[ or, 
more recentely, in Teyssandier (22j . Moreover, the ana- 
lytical solution for light propagation in the gravitational 
field of one spherically symmetric body in the framework 
of post-post-Newtonian (ppN) development has been re- 
analyzed [231 ]. Analytical formulas have been compared 
with high-accuracy numerical integrations of the geodetic 
equations, the latter based on an approximation of the 
exact analytical solution. The authors found a deviation 
between the standard post-Newtonian approach and the 
high-accuracy numerical solution of the geodetic equa- 
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tions which amounts to 16 //-arcsecond. As reported by 
Klioner and Zschocke [23] ]. detailed analysis has shown 
that the error is of post-post-Newtonian order in contra- 
diction with the analytical estimate. To clarify the prob- 
lem, standard post-Newtonian and post-post-Newtonian 
solutions have been derived by retaining terms of rele- 
vant analytical orders of magnitude. For each individual 
term in the relevant formulas exact analytical upper es- 
timates have been found and a compact analytical solu- 
tion for light propagation is derived up to 1 /i-arcsecond 
accuracy. The enhanced post-post-Newtonian terms cal- 
culated by Klioner and Zschocke, are confirmed recently 
by Teyssandier [24[ 

On the other hand, Kopeikin and Schafer 
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the pM approximation, solve Einstein's equations in the 
linear regime by expressing the perturbed part of the 
metric tensor in terms of retarded Licnard-Weichert po- 
tentials. Later, Kopeikin and Mashhoon [13J included all 
the relativistic effects related to the gravitomagnetic field 
produced by the translational velocity/spin-dependent 
metric terms. In these works a special techni que of in- 
tegration of the equation of light propagation [25| with 
retarded time argument is extended both outside and in- 
side of a gravitating system of massive point-like particles 
moving along arbitrary world lines. The null geodesic 
is rewritten as function of two independent parameters, 
one related to the Euclidean scalar product between the 
light ray direction and its unperturbed trajectory, and 
the other derived as the constant impact parameter (26j 
of the unperturbed trajectory of the light ray. Then, the 
light trajectory is again traced as a straight line plus in- 
tegrals, each containing the perturbations encountered, 
from a source at an arbitrary distance from an observer 
located within the n-body system. 

As a consequence, one more difference between the 
pN/pM approximations appears in the computation of 
the light deflection contributions: in the pN scheme this 
is done by combining the external and internal solution 
with the technique of asymptotic matching, while the pM 
method utilizes in addition a semi-analytical integration 
of the equation of light propagation from the observer to 
the source with respect to the retarded time considered 
as argument in the metric tensor. The relativistic per- 
turbations remain expressed in terms of integrals (i.e. a 
semi-analytical solution) which can be solved analytically 
in specific cases (i.e. not expressed as a general solution 
of the differential equation, but as particular solutions 
for some effects on the light propagation) once the mo- 
tion of the gravitating bodies is prescribed. The pM ap- 
proximation is also used by Teyssandier and Le Poncin- 
Lafitte [13] that found a way to bypass the solution of 
the null geodesic in form of differential equations. In this 
work the authors elaborate the light trajectory with the 
powerful tool of the Time Transfer Function [28[ , worked 
out throughout a recursive procedure for expanding the 
Synge's World Function, a covariant method that is con- 
ceptually much closer to that of RAMOD, the approach 
discussed in the present work. 



RAMOD stands originally for Relativistic Astrometric 
MODel, conceived to solve the inverse ray-tracing prob- 
lem in a general relativistic framework not constrained 
by a priori approximations. RAMOD is, actually, a fam- 
ily of models of increasing intrinsic accuracy all based on 
the geometry of curved manifolds ( [29|, H3| and reference 
there-in) where, as this paper proves, light propagation 
can be expressed in a general relativistic context. An 
analytical solution of the null geodesic exists in the case 
of the Parametrized-Post-Newtonian Schwarzschild met- 
ric 0, HH, which, just considering the spherical mass of 
the Sun, can be implemented for the Gaia mission (Esa, 
Turon et al. (32[) and could prove the dilaton-runaway 
scenario (33|. At present, the RAMOD full solution re- 
quires the numerical integration of a set of coupled non 
linear differential equations, called "master equations" , 
for tracing back the light trajectory to the initial position 
of the star, and which naturally entangles the contribu- 
tions of the curvature of the background geometry. A 
solution of this system of differential equations contains 
all the relativistic perturbations suffered by the photon 
along its trajectory due to the intervening gravitational 
fields. The boundary condition are usually fixed by the 
required physical measurement (34j . In this context, a 
semi-analytical solution was found by de Felice and Preti 
35]. 

This work shows how RAMOD can naturally recover 
the seminal work of Kopeikin and Schafer [12] and 
Klioner [l4| ■ While the third one stems from the gen- 
eral formulation of the second by conforming to its ba- 
sic conception [l8j], RAMOD responds to an inherently 
different strategy. As far as the invers light ray trac- 
ing problem is concerned, in the Kopeikin and Shaicr 
and Klioner approaches the light ray is reconstructed as 
a sum of terms which allows for a direct evaluation of 
the individual relativistic effects, induced by the grav- 
itating bodies, suffered by the light on its way to the 
observer. RAMOD, instead, aims to determine a full so- 
lution for the light trajectory which naturally include, in 
a curved space-time, all the individual effects; the latters 
are somewhat hidden in the covariant formalism of this 
approach and directly contribute to the solutions of the 
master equations. Evidently, any specific effect one is in- 
terested to explore can be independently deduced from 
our formalism as a branch output, once we adapt the 
model to a required specific case. A proof of this is given 
in Crosta and Vecchiato (36| . where a first comparison 
between RAMOD and GREM (Gaia RElativistic Model, 
jl4j, which employes the pN formulation) was carried out 
via the extrapolation of the aberrational term in the " lo- 
cal" light direction, i.e., at the observer. 

Now, if we consider a gravitational field generated by 
a n-body system, the required accuracy in back trac- 
ing the light ray should be consistent with the geome- 
try induced by the system itself. The accuracy is gauged 
by the virial theorem, i.e., by the ratio v/c, where v is 
the typical velocity of a body of the system. In this re- 
spect, RAMOD3 [29| is developed up to the order (v/c) 2 , 
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namely, it is based on a static geometry of the space-time, 
whereas the most recent RAMOD4 [30jj is the extension 
of the predecessor to include a dynamical geometry (i.e., 
all metric terms up to the order (v/c) 3 ), which is ac- 
complished by including also the vorticity contribution 
proportional to the off-diagonal term goi of the metric. 
Both models introduce the retarded contributions to the 
metric due to the gravitational influence of the sources 
on the light path, according to the natural space-time 
structure. In section IIIIl we discuss how to build an 
appropriate geometrical set-up in order to handle light 
propagation through such a space-time. In section HVl we 
present the building process of the differential equations, 
i.e., the master equations, from which the light trajec- 
tory from an observer, located within the n-body system, 
to a distant star can be reconstructed. We distinguish 
two cases: (i) the static space-time (RAMOD3), and (ii) 
its dynamical extension (RAMOD4). An application of 
RAMOD3 to the euclidean metric is presented in section 
fVl while section IVII shows how this last computation al- 
lows to parametrize the RAMOD3 master equation as 
done in [12[ ■ Finally, section I VIII shows how the veloc- 
ity of the bodies can enter the metric term simply as a 
"retarted distance" corrections in RAMOD3-like models. 



II. NOTATIONS 

In this paper the following notations will be used: 

• components of vectorial quantities are indicated 
with indexes (no bold symbols), where the latin in- 
dex stands for 1,2,3 and the greek ones for 0,1,2,3; 

• regular bold indicates four- vector (e.g. u); 

• indexes are raised and lowered with the metric g a p ; 

• rfrii or n l n l stands for the scalar product with re- 
spect to the euclidean metric 5ij, whereas l a l a with 
respect the metric g a p; 

• c indicates the fundamental speed (equal to the 
speed of light in vacuum); 

• A quantity like T a , represents the tangent vector 
to a curve T; 

• P(u) Q( 3 represents the operator which projects or- 
thogonally to its argument. 

• x' a or x a indicate generic four-coordinates, whereas 
£ a the Lie-transported ones; 

• dot applied over coordinates means derivative with 
respect to the coordinate time; 

• partial derivative with respect to the coordinates 
are usually indicated by a comma. Instead, the 
symbol d is adopted in cases where the formulas 
need to be more explicit and avoid confusion; co- 
variant derivative are indicated with V. 



III. SETTING UP THE GEOMETRY 

The RAMOD framework is based on the small cur- 
vature limit, for which the curvature of the background 
geometry is sufficiently small to neglect non linear terms. 
The degree of this approximation has in turn to be spe- 
cialized to the particular case one wishes to model. The 
first step is to identify the gravitational sources and then 
fix the background geometry. 

Let the space-time be generated by a weakly relativis- 
ts gravitationally bound system; in this case the metric 
is 

g a p = Vap + h a/ 3 +0(h 2 ), (1) 

where rj a p is the flat Minkowskian metric and the h a p's 
describe effects generated by the bodies of the system and 
are small in the sense that \h a p\ 1. From the virial 
theorem all forms of energy density within the system 
must not exceed the maximum amount of the gravita- 
tional potential in it, say, U. So, the energy balance 
requires that \h a p\ < U/c 2 ~ v 2 /c 2 , where v is the av- 
erage relative velocity within the system [37j . Since the 
latter is weakly relativistic [38f . the h a p's are at least of 
the order of (v/c) 2 and the level of accuracy, to which is 
expected to extend the calculations, is fixed by the order 
of the small quantity e ~ {v/c). In summary, the pertur- 
bation tensor h a p contributes with terms even in e to poo 
and gij (lowest order e 2 ) and with terms odd in e to goi 
(lowest order e 3 , Misner et al. [3lh ; its spatial variations 
is on the order of |ft. Q /3|, while its time variation is on 
the order of e|/i a| a|. Clearly the metric form ([1) is pre- 
served under to infinitesimal coordinate transformations 
of order \h\. 

Under these hypothesis, in order to describe the space- 
time evolution of the system, let us introduce a family of 
physical observers, i.e. a time-like congruence of curves, 
and consider its vorticity, which measures how a world- 
line of an observer rotates around a neighboring one. If 
an open set of the space-time manifold admits a vorticity- 
free congruence of lines, then it can be foliated (Frobenius 
theorem, [39H40j). 

Given a general coordinate space-time representation 
(x l , f')i=i,2,3, the associated metric being g' aj3 (x' ,t'), 
there exists a family of three-dimensional hypersurfaccs 
S(r), or slice, described by the equation: 

t(x' ,t') = constant, (2) 

where t(x' ,t') is a real, smooth and differentiable func- 
tion of the coordinates. A space- like foliation implies that 
a unitary one-form u' a exists which is everywhere propor- 
tional to the gradient of r (figure [T]), namely (3ll. l39l| : 

"a = (3) 

where 

g'^u' a u's = -1, (4) 
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and 



e v/ = [-</ a/ wx^<)r 1/2 . 
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FIG. 1. A coordinate grid (x',t') and the hypersurface 
r(x',i') = constant. 

Evidently the curves tangent to the vector field, i.e. 
u a = g' a Pu'p, form a congruence of time-like lines C u ', 
everywhere orthogonal to the slices Sir) . Let us compute 
the vorticity as the tensor quantity [39j 

uj af) = P' a p P' fj a u' [pij] , (6) 

where square brackets mean anti-symmetrization and 

P'«p{u , )=g' a p + u' a u' fj (7) 

is a projector operator on S(t). Since Uj Q ~, results in 
u'^d^tp' , we find that the congruence C u > is vorticity- 
free only up to (v/c) 2 order; namely, 



(8) 



Then, we may affirm that considering only terms up to 
e 2 r~j (v/c) 2 , it is the same as assuming that structure 
of space-time is static. So, does the one-form u' a have 
a physical meaning? Let us introduce a new coordinate 
system (x°,x l ) such that: 



where if> = ip(x°,x) is a unitary function (u a u a = — 1), 
a is the parameter on the normals, and 



u a (x°,x l ) = -^—u (x',t') 



dr 
"dx^ 



dx° 

>'(x',t') _ _ S O e ^(x°y) _ 



Clearly it is 



da 
dr~ 



(-«") 



-1/2 



(11) 



(12) 



which means that the parameter a runs uniformly with 
the coordinate time r. The congruence C u , however, 
does not Lie-transport the spatial coordinates x % . From 
© and (|A2[) . the line element can be cast into the form: 

ds 2 = -{Ndx°) 2 +g l3 (dx l - N i dx°)(dx j ~ N j dx°) (13) 

where A" = ±l/u° and AP — u l /u° are termed respec- 
tively lapse and s/ii/£ functions (or factors) [3l|. In fact, 
choosing at time n, say, an arbitrary event T 3 , labeled 
by the coordinates (ti , a^) , the unique normal through 
that point will intersect the slice S(ti + At) at a point 
V' with spatial coordinates which are shifted from the 
initial ones by the amount (see appendix A): 



Ax 1 = 



At 



N i (T')dT'. 



(14) 



The spatial coordinate transformation to be applied, in 
order to make the shift factor vanish and assure that the 
congruence of normals will now Lie-transport the spatial 
coordinates, is: 



d£ = dx l - N l dx° 
d£° = dx° = dr, 



in which we impose: 



fjpa 



(15) 



(16) 



Under the above transformation, we clearly have for a 
0: 



fV + f^Mo ■ = -.-■ 

OT OX 



x l = x l (x' ,t') 



(9) 



where x 1 are the space-like coordinates defined on each 
slice. From (j9|), the congruence of normals now have the 
following tangent vector field: 



'u a (x°,x l ) = -S° a e^ 
u a (x°,x % ) 



dx a 
da 



(10) 



= ^9 



no 



from (O; and for a — i 



1T U + 1T1 U = TT9 + 7T19 = °- 

OT OX 1 OT 0X3 



But it is: 



r 



dx a - 
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which implies 



-00 



(17) 



90i 



0. 



(18) 



Then the vector field tangent to the congruence Cq 
which transports the spatial coordinate is: 

(19) 



da 

or 



Let us summarize. First, the vorticity is proportional 
to the goi term of the metric, the lowest order established 
by equation [5] second, we must bear in mind that, even 
if metric form ([l} holds up to infinitesimal coordinate 
transformations, the condition that the shift factor is zero 
is preserved under the above transformation only to the 
order e 2 . In fact, to the order of e 3 , the terms g' Qi will not 
remain zero if they were so initially. We may affirm that 
considering only terms up to e 2 ~ (v/c) 2 , it is the same 
as assuming that structure of space-time is static, since a 
static space-time is a stationary spacetime in which the 
timelike Killing vector field has vanishing vorticity, or 
equivalently (by the Frobenius theorem) is hypersurface 
orthogonal, i.e. admits a family of spacelike surfaces of 
constant time ( [HI, |39| ) . The parameter a along the 
Killing congruence Cq, such that «° = d£ a /dcr, is the 
proper time of the physical observers who Lie-transport 
the spatial coordinates. Any hypersurfaces, at each dif- 
ferent coordinate time r, can be considered the rest space 
of the observer u. One can require that the world line of 
the center of mass of the body's system, chosen as ori- 
gin of coordinates of the gravitational bounded system, 
belongs to this congruence while the world lines of the 
bodies would differ from the curves of the congruence by 
an amount which depends on the local spatial velocity 
relative to the center of mass; however, at the e 2 order, 
i.e. neglecting the relative motion of the gravity sources, 
one can assume that also their world lines belong to that 
congruence and the geometry that each photon feels is, 
then, identified with the metric g a /3- To the order of e 3 
the shift factor should be retained and so also terms pro- 
portional to hoi (eq. |S]) . Since we require that the back- 
ground geometry is weakly relativistic, we can assume 
that the vorticity is free up to the e 3 order only locally; 
then the space-time still admits a family of hypersurfaces 
r ^constant with normals forming a congruence of curves 
with tangent field u given by (|10l) . The given coordinate 
system (x l , r), where x l are fixed on each slice, is still 
centered to the baricenter of the n-body system, but, on 
the other hand, each slice is not the local rest-space of 
the observer u. These normals do not in general Lie- 
transport the spatial coordinates x l ; instead they vary 
along the normals according to the shift law (fT4)l . Any 
of these observers can be considered at rest with respect 
to the coordinates x l only locally, and for this reason u 
is called local barycentric observer. 



IV. FROM THE NULL GEODESIC TO THE 
MASTER EQUATIONS 

In the following sections we drop the tilde for the met- 
ric. A photon travelling from a distant star to an ob- 
server, located within the bounded system, would see 
the space-time as a time development of the r = con- 
stant slices. Let us now consider a null geodesic Tfc with 
tangent vector field k a = d£ a /dX which satisfies the fol- 
lowing equations: 



k k(y^ - — 

rlh a 

^ + r^ = 0; 



(20) 
(21) 



here A is a real parameter on Tfc and r" CT are the con- 
nection coefficients of the given metric. Assume that the 
trajectory starts at a point T 3 * on a slice 5(r*) (say) and 
with spatial coordinates £*. The null geodesic crosses 
each slice S(r) at a point with coordinates = £ 1 (A(t)); 
but this point also belongs to the unique normal to the 
slice S(t), crossing it with a value of the parameter 
a = ct(£*(A), r) = cr^i( A )(r) (figure^. Let us now de- 




FIG. 2. The point £*(A(t)) belongs to the null geodesic but 
also to the unique normal to the slice S(t) crossing it with a 
value of the parameter c^imC 1 ")- 

fine the one-parameter local diffeomorphism (39| : 



»A<T 



,(t)) 



T t nS(T)^S(r„) (22) 



which maps each point on the null geodesic Tfc with the 
point on the slice S(t ) which one gets to by moving along 
the unique normal through the point Tfc (A) n S(r), by a 
parameter distance Act = C£»(A)( T o) — (T i i (X)(' r ) (figure[3]). 
Since the spatial coordinates £ l are Lie-transported along 
the normals to the slices, then the points in S(tq), which 
are images of those on the null geodesic under <^a<t, have 
coordinates (<M<x(Yfc((A) n S(t))) 1 = £\ The curve in 
S(tq) which is the image of Tfc under 0a<t, is: 



>(-)) 



o Tfc 



(23) 
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with tangent vector |3£ 



This coincides with the projection operation on the 
rest-space of the observer u in any point of the mapped 
trajectory on the silde: 



£ a = P(u)%k p , 



(25) 



hence the curve T is the spatial projection of the null 
geodesic on the slice S{t ) at the time of observation; 
this curve will be denoted as Ti and is naturally param- 
eterized by A. Then from equation (j2"4")l by using (TT9"]) 
(where u a u a = —1) or (f2"5j) . it follows that: 



£ a = k a 



(26) 



Clearly by setting a = and from equation (|T9"f we de- 
duce £° = 0; moreover, l a is space-like as expected as it 
is a projection on slice S(tq) and lies everywhere in it. 
Since each point of T ' g is the image under the diffeomor- 
phism </>A(t of a point Tfc D S(t), it is more convenient 
to label the points of Ti with the value of the parameter 
°f i (A)( r ) which, as we have already said, uniquely iden- 
tifies that point on the normal to the slice S(t) which 
contains it. Hence, being 



da = —(u a k a )d\ , 
we define the new tangent vector field: 



In the same way we denote 



(u^kp) 



(27) 



(28) 



k a 



so that 



which implies: 



k a = l a + u° 



i a i n = 1. 



(30) 
(31) 



Note that the parameter 0£<(A) is also the one which 
makes the tangent to the curve T unitary. To short the 



notations, in what follows we denote a 



as a. We can 



now write the differential equation which is satisfied by 
the vector field t. From {27]) and equations (|2l"j) 

becomes (see appendix [B| : 



da 



du a 

i 

da 



a + u Q )(^{2 r v r u,3 + efirvrup) 

+ = 0. (32) 






u 








Ti 

















FIG. 3. The curve T is the image of the null geodesic Tfc 
under the diffeomorfism 4>& a . 

In this equation the quantity £^E T V ' T Up can be written 
explicitly in terms of the expansion O ptT of Cq [29| . as: 

it + ^~ ~ ^ + u a )(^u fj + e p J»n 

da da 

+ T a M {F + u'^iP + u 1 ) = (33) 

where 9 pCT = P°PgV { a up) |39| . Since the only non van- 
ishing components of the expansion are Qij — (l/2)dohij, 
the expansion vanishes identically as consequence of the 
assumption to neglect time variations of the metric. 
From this condition equation, after some algebra, (|33[) 
becomes to the required order (see appendix [Cjl : 



A pa Jf.Q 

==- + ^- - {£ a + u a )(i%)+ rs 7 (P + + iP) 

da da 



(29) If a — 0, equation (j3"2")l leads to 



dF 

da 



= 0. 
(34) 

assuring that 



condition £° — holds true all along the curve T; if a = k 
equation (|32|) gives the set of differential equations that 
we need to integrate to identify the star or the physical ef- 
fects related to light propagation (see details in appendix 
El gl| and H3): 



dt k 
da 



1 



1 



-r [r/l'-uj - S ks \h aJii - -h lhS ) i'i ' 



(35) 



These equations determine light propagation in the 
static case, here and after called RAMOD3 master equa- 
tions. Let us remark, that this definition does not stand 
for the "master equation" used in classical or quantum 
physics; in this context it represents a set of first-order 
nonlinear differential equation describing evolution of the 
spatial light direction components according to the pre- 
scribed geometry. In the dynamical case, any "local" 
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observers u which are at rest with respect to the coor- 
dinates x % locally, once intersected by the null ray will 
see the light signal along a spatial direction I in his rest 
space given by (more details can be found in de Felice 
et al. 13011) 



£ a = P^ Q (ii)fc' 3 (i 



(36) 



where P(u) a p — 5^ + u a U[j this time is the operator 
which projects orthogonally to the u's. Following the 
same reasoning of the static case, we parameterize the 
curve T with the parameter er which makes unitary the 
locally projected vector field I which we term again t, 
so l a l a = 1. The differential equation of the null ray, 
written in terms of £ a , is still given by (|3"2"j) . the only 
difference being that hoi ^ 0, /ioo,o 7^ 0, and hijja 0. 
After some algebra, equation (|32|) now reads: 



Assume that the mapped spatial photon trajectory T^, 
in the slice t — tq, belongs to a normal neighbourhood 
of CM. By definition then, there exist a unique geodesic 
(in any case spacelike geodesic) T' connecting CM to 
any point V of Tg. Denoting with <t(t, £*) the parameter 
along T^, with A the parameter along each of the curves 
T stemming from CM and such that Aqm = 0, and with 
<p the homeomorphism which defines the chart containing 
V in the normal neighbourhood of CM, we have [3!| : 

P = T(a) = T'(A) = ^ 1 (f(a)). 

The coordinate can be expressed as: 

C = ^) = b°T(<r)] i = [^T'(A)] i . 

Applying to both members of the above equation the map 
T _1 o (p^ 1 , we obtain : 



da 



dP k 1 - - - 



(37) 
(38) 



+ ^thoo.i + P {h k o,i + hki,o - hoi,k) 



1 



hook — £ Phoifi + hka.o — 0. 
z 

These equations are named RAMOD4 master equations. 
Note that in RAMOD4, with respect to RAMOD3, we 
have one more differential equation for the £° component. 
Term like /ifco.o is redundant at the order of e 3 , but we 
keep it in order to compare the geodesic equations, as it 
will be clear in section 6. 



PARAMETRIZED MAPPED 
TRAJECTORIES 



A=(T'- 1 o^ 1 )o( v ,oT(a))=A(f(a)), 
namely the A implicit dependence on a at the point V 
with coordinates £ l . Denoting as — 



dA 



the unit tan- 



gent to the curve T , we define the length of T between 
CM and V as the quantity: 



L = 



/ (CC0 1/2 rfA'^A(a). 
Jo 



(39) 



The point of closest approach of the photon to the center 
of mass, will correspond to the minimum of L. Thus, 
differentiating (|39|) with respect to a, we have 



and at the point of closest approach with parameter a 

(d i% = o. (40) 



Let us fix the origin of the coordinates and consider it 
as the center of mass (CM) of the matter distribution of 
the n-body system. Let us, also, choice that its spatial 
coordinates belong to the congruence of curves Cq. 




Let us set A(<r) = £ as the impact parameter which rep- 
resents the distance from V(a) to CM; this is a constant 
which labels each mapped photon trajectory on any given 
slice t = tq. Recalling that £ I (A(cr)) = £ l (cr) at any point 
V of intersection of the curves T with the curve T , then, 
at CM, where Ccm = and Acm = 0, we have to the 
first order in A: 



OA \ ' ^ ( ^ 



6 A 



(Co 



8 A 



(41) 



(CM->P) 

where (£§) (cm->v) re ^ ers to tr±e tangent vector at the 
origin of the "unique" spacelike geodesic connecting CM, 
the origin, to the point V . From (|41 [) we clearly have: 



FIG. 4. Mapped photon trajectory in a slice to and normal 

neighbourhood of the origin of the spatial coordinates £ l (CM) C(A) = (Co) 



dA = {<*) 



A (42) 
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Let us again introduce a new parameter on T, namely 



then, at a point V on T with parameter a + Sf 



But from (j42|), at A(a), it is also 



e (A(a)) = (a) 



hence, at any point V on T, we obtain 



= + 



frdf 



(43) 



(44) 



(45) 



i.e., the coordinate path of the photon on the slice To as 
function of the two parameter £ l and f. Formula (j45|) has 
been obtained for a mapped trajectory and it contains the 
local spatial direction in integral form. In this sense gen- 
eralizes the parameterization of Kopcikin and Mashhoon 
[l3| (and any references therein) , where, instead, the Eu- 
clidean scalar product is used (i.e. tli — SijftP). Hence, 
if the euclidean scalar product is applied and in case of 
constant light direction, equation (|4"5j) is identically the 
same. However, by calculating the square modulus of 
(1431) we get 



f 2 . 



(46) 



The reader should bear in mind that equation (1431) has 
its validity only up to e 2 and can be traced onto the 
hypersurface where the photon's trajectory is mapped. 



VI. NULL GEODESICS IN COMPARISON 



The quantity £ a (eq. [23) in RAMOD is the unitary 
four- vector representing the local line-of-sight of the pho- 
ton as measured by the local observer u, at rest with re- 
spect to the barycenter of the n-body system. In short, 
it represents a physical entity. Scope of this section is 
to demonstrate, first, that substituting its coordinate 
counterpart, within the appropriate approximations, to 
equation (13"5|) of RAMOD4 is equivalent to recover light 
geodesic equation in the first post-Minkowskian regime 
adopted in Kop eikin and Mashhoon l_3j or Kopeikin 
and Schafer jT2[- Second, the above authors parametrize 
the coordinate expression of the null geodesic to solve 
the light propagation problem in many physical context, 
as the literature largely shows, but, as far as RAMOD 
is concerned, only equation (|3"5"|) of RAMOD3 recovers 
the parametrized one (37) in Kopeikin and Mashhoon 
13] (equivalent to equation (19) in Kopeikin and Schafer 
12]). Nevertheless, let us use the definition of £ a ; since 



it is 



k' 



u a k a ' u Q k°[-l + h 00 + h 0l (k i /k )Y 



with u° = (— goo)^ 1 ^ 2 an d k l /k° — dx t /dx°, the spatial 
coordinate components of £ % result (in what follows we 
assume c = 1): 



1 



i i =x i (-g 00 ) 1/2 (l-h 00 -h 0i x i y 

2 ho ° 



h 0l x' 



Oh 2 



(47) 



Let us transform equation (|3"51) of RAMOD4 into its coor 
dinate expression and consider only its spatial part. Ne- 
glecting all the contributions non linear in h, from (|47l 
and (1121) the left-hand side results: 



dl k 
do 



d 



1 



x k {\ + -h 00 + h 0i x l ) 



1 



x K (l + -h 00 



1 dh 



2 dr 



- hoiX 1 
dh 0i . 

~dV x 



= x k + ± k 



hoix* 



2^oo,o+ 



hoiji'x 3 +h 0h0 x l ) +0(h 2 ), 
while the right-hand-side transforms as: 



(48) 



x x^ h^jQ x x? \^hkj,i 
^x k x l h 00 , l - x l (hko,i + hkifi 
r/ioo,fe + x k x % hoi^ — hko,o 



) + 



(49) 



Equating expressions (|48| and (|49|) we get the approx- 
imated null geodesic equation (40) in Kopeikin and 
Schafer (l2l ] (and references therein), namely in coordi- 
nate form: 



1 



hoo.k — h()k,0 — -hoonx — hkinX 1 



2 

{hok.i 



hoi,k)x' - h 00tl x k x l 



-hij. k ) x l x J + 



(50) 



We remind that the term x k is of the order of h because 
of the order of the geodesic equation itself. Actually, 
this last result should be expected, since both equations, 
(P]l of RAMOD4 and (40) in Kopeikin and Schafer [13], 
are deduced from the null geodesic (|2"Tj) in a weakly field 
regime. Once obtained this equivalence at the h order 
by using the coordinate expression (|47[) of £, in principle 
one could adopt the same parameterization done, for in- 
stance, by Kopeikin and Mashhoon [13| (equation (36), 
or Q, equation (7) and reference therein) in order to 
find a solution of the master equation in the RAMOD 
framework applicable to all of the physical cases already 
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presented in the literature in this context. These parame- 
ters are: (i) r, defined as the Euclidean scalar product be- 
tween the light ray vector and its unperturbed trajectory 
(eq. (13) in [13j , for example), treated as a non-afhne pa- 
rameter, and (ii) £ = the constant impact parameter 
of the unperturbed trajectory of the light ray. But as sec- 
tion [V] shows, this kind of parameterization in RAMOD 
is possible only in a static space-time where null geodesic 
can be entirely mapped into a hypersurface simultaneous 
to the time of observation, i.e. at the e 2 level of accuracy 
in the linear regime. Then, consistently to the previous 
reasoning, we need to check if RAMOD3 master equa- 
tion can be transformed into eq. (36) of Kopeikin and 
Mashhoon Therefore, let us adopt the parameteri- 
zation (l45j) in case of a constant light direction, say £ ; 
moreover, let us introduce the impact parameter p as 
done in |25j |. i.e. assume the equivalence of the two pa- 
rameterizations as a consequence of the following change 
of coordinates: 



i k = i'!c 

? = r. 



£ k r + E k , 



(51) 



In this transformation the impact parameter p is as- 
sociated to the coordinate projected orthogonally (with 
respect to the light direction) by the operator PJ" = 

S'p— l a lp; E k represents a correction to the coordinates of 
the same order of the perturbation term \h\ of the metric, 
and df = da. Once adopted such a new coordinate sys- 
tem, the partial derivatives should also change according 
to the following rule: 



d _ dp d d£° d 
W ~ W~dp~ + WW 



^ = ^^7 + ^^ = -P/^+3V (52) 



d_ 



The next step is to compute the new components of the 
metric h a p according to transformations (|5ip : 



dp dp 



~ Q^O g£0 h ™ f Q£0 d £0 "PO n g ^ g ^.„ 

= h 00 -2e h 0p +e £fh gp + o{h 2 ) , 



(53) 



_ dp dp t dp dp t dp dp i 

hoi WW 00 + WW p0 + WW q + 

ww hpq = p * hoq ~ ilP?hpq + ° {h2) ' (54) 



and 



dpdp t dp dp* dpdp~ 

hij - WW + WW p0 + WW q + 



dp dp 

WW 



h pq = Pfpfh pq + 0{h 2 ). 



(55) 



We can recast the RAMOD3 master equation by applyng 
formula (|47l) and keeping in mind that, all along the 
mapped photon trajectory, we are not allowed to neglect 



the derivative of the metric with respect to a, or the new 
parameter f, since they vary, for each Lie-transported 
spatial coordinate, according the one-parameter local dif- 
feomorphism (f2"3"f 

~ Vfc - \k 00 ,4 k ? - {hki,j - CP, (56) 

which is equivalent, setting £q as a constant, to: 



d 2 E k 1^„, ldftoo, fc 1 



df 2 2 



~Pk h 00,q 



100 ok x p«i, ok pi 
2 UT 2 



[Pjhu, q - \p q k h lhq ^j W (57) 

Now, considering also that dhoo/df — P^hoo^£ + ftoo,f > 
we can eliminate the projected term along the corre- 
sponding coordinate direction and, after substituting the 
new metric terms, find: 



d 2 E k 1 



df 2 2 
1 
2 



(h Q0 - 2£ p h p0 + £ p £fh p 
\ (h 00 - 2£ P Q h P o + £ P £fh Pq ) 



l pPp1u \ fifl 
2 i j PI j „ ^0*0' 



d 2 E k 1 



df 2 



- (hoo - 2£ p h p0 + £ p £fh 
^(hoo-2£ p h p0 +£ p £ q h pq ^£ h . (58) 



The last equation is not yet comparable to the one in 
Kopeikin and Mashhoon [13|, but we have to consider 
the fact that the non-diagonal term like goi in the static 
case are null, so from (|5l)l 



(59) 



Pfho P — £\P p h pq 



- t p £i>ho P = -hi + £lk q - £ q o£ P £ 0l h pq . (60) 

Once replaced the right-hand-side of eq. (|6"0"|) in eq. (|58l) , 
we straightly get: 



d 2 E k If- p ~ p ~ 

-TT2- ~ 2 \ h o° - 2£%h p o + £ P o£ q oh mj k 

\hwA ~ h ok + £ q h kq - ^£ k o£ p £ q h pq 



(61) 



which is exactly the same equation (37) Kopeikin and 
Mashhoon 13]. Actually, at this stage, these authors 
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introduce the four-dimensional isotropic vector k a 
( — 1, k l ) in order to get the following equation: 



d 2 E k ^ 
df 2 ~ 

— k k^ h a ? 

2 ap ,/c 



(62) 



-h 00 k k 



k hfcfy ~*2i^ k^ k^ hp 



In conclusion, as far as RAMOD framework is concerned, 
equation (|61[) can be considered as an ordinary, second 
order differential equation for the perturbation 5 fc in vari- 
ables valid in the domain where e 2 accuracy holds, i.e. 
where RAMOD3-like master equations must be used. 



VII. BODY'S VELOCITY TERMS IN THE 
METRIC 

The integration of the master equations requires to cal- 
culate the metric coefficients h a p, which depends on to 



the retarded distance rt a \ as discussed in 



and [3C 



This means that we have to compute the spatial coordi- 
nate distance ri a \ from the points on the photon trajec- 
tory to the barycenter of the a-th gravity source at the 
appropriate retarded time and up to the required accu- 
racy. Exploiting the mapping procedure we have used 
to identify the spatial photon trajectory on the slice at 
the observation time, we can map on the same slice the 
trajectory of each body of the system. We shall term 
£*(cr) the spatial coordinates along the spatial photon 
path X and x % (a) those on the spatial path of the a-th 
source, a being the parameter along its spacetime tra- 
jectory (see figure [5]) . Evidently, the metric coefficients 
at each point of the photon's spatial path with param- 
eter tr(r) corresponding to the photon's spatial position 
at the coordinate time r, are determined by each source 
of gravity when it was located at a point on its spatial 
path with parameter ct(t'), where r' is the retarded time: 
t 1 = t — r (c = 1). The coordinate distances is given by 
(we drop the suffix (a)): 

r = \e{a(r))-x\a{r'))\ (1 = 1,2,3). (63) 

Since all the functions here are smooth and differentiablc, 
we can expand the coordinates of the gravitating bodies 
in Taylor series around their position at time r as: 

x\d{T + 5T)) =x 1 {o{t))+ Sr + --- (64) 



As reported in papers [29j, |30j , the right-hand-side of (|63|) 
is equivalent then to (c — 1): 

e(a(r)) - x\a{r')) = f(a(r)) - x\a{r)) + (65) 



where v a is the spatial velocity of the source in the rest 
frame of the barycenter. If we want our model be ac- 
curate to e 3 , then it suffices that the retarded distance 



r contributed to the gravitational potentials, which we 
recall are at the lowest order e 2 , at least to the order of 
e. Hence equation ([63]) leads to: 



v j d( + 0[ 



where 



and 



ro 



Ae = e(v(T))-x\a(T)). 



(66) 



(67) 



(68) 



Evidently, to the order of e 2 (static geometry) the contri- 
bution by the relative velocity of the gravitating sources 
can be neglected. Indeed, in the static case one can fur- 
ther expand the retarded distance in order to keep the 
terms depending on the planet's velocity up to the de- 
sired accuracy. Obviously the gravitational field does not 
vary as photon moves from the star to the observer, since 
terms goi are null and the time derivatives of the metric 
are of the order of e 3 , therefore, the effects due body's 
velocity cannot be related to a dynamical change of the 
space-time. Actually, the positions of the bodies at dif- 
ferent values of the parameter ct(t) can be recorded as 
subsequent snapshots into the mapped trajectories and 
deduced as postponed corrections in the reconstruction 
of the photon's path. 



VIII. CONCLUSIONS 

Most of the information, from the Universe around us, 
comes in form of electromagnetic waves. While propa- 
gating, light interacts with the curvature generated by 
the encountered mass. Then, it is of utmost importance 
to reconstruct the whole light path precisely as much as 
the varying curvature. For such a reason, the covariant 
RAMOD framework aims to use and preserve as much as 
possible all the physical quantities needed to reconstruct 
the history of the photon. As a matter of fact, modelling 
light propagation is intrinsically connected to the identi- 
fication of the geometry where photons naturally move. 
This process can substantially influence the solution of 
such a problem, since, according to the physical environ- 
ment where observations take place, appropriate metric 
terms should be retained or ruled out in the affine con- 
nection that contributes to make a null geodesic explicit 
and integrable. 

The result reported here, in the case of RAMOD-like 
models, represents an improvement in the understanding 
of how to handle the null geodesic within the accuracy 
requirements due to the geometry, i.e. the physics of a 
n-body bound system. To avoid confusion, the reader 
should bear in mind that the present work does not deal 
with the description of the n-body dynamics; it consid- 
ers, instead, how the set-up of the space-time for light 
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FIG. 5. Mapped curves related to the planet's retarded po- 
sitions and retarded distance r. The pink line represents the 
mapped trajectory of the photon, here represented at time r 
with a light cone. Also the coordinates of the planet x l (a(r)) 
can be Lie-transported along the congruence of curves and 
mapped on the same slice (blue line). The contribution to 
the metric coefficient at each point of the photon's spatial 
position at coordinate time r is determined by the source of 
gravity when it was located at a point on its spatial path, 
along the light cone, at retarded time r'; r, in this picture, 
represents the spatial coordinate distance on slice S (to) be- 
tween the mapped photon's position at coordinates time r 
and the mapped source's location at coordinate time r'. 



propagation and detection changes without or with the 
vorticity term, being the latter related to the congruence 
of curves that foliates the space-time and allows to de- 
fine a Lie-transported coordinate system. The vorticity 
term cannot be neglected at the the order of e 3 : the pos- 
sibility to ignore it locally can be applied only to a small 
neighborhood with respect to the scale of the vorticity 
itself. In other words, we have to scrutinize if a measure- 
ment can be considered local or not with respect to the 
curvature [ioj . 

Within the scale of the Solar System, then, there is 
none slice which extends from the observer up to the star 
that emits light. This justifies why only RAMOD3 re- 
covers the paramctrization used in Kopeikin and Schafcr 
[l2l |. namley when the geometry allows to define a rest- 
space of an observer everywhere. In this regard, this 
paper proves that in RAMOD3, i.e., a static space-time 
model, the parameterization (I45[) is appropriate to solve 
the light ray tracing problem. Formula (|4"5j) has its valid- 
ity in e 2 regime for a mapped trajectory and it contains 
the local spatial direction in integral form. It general- 
izes the parameterization used in Kopeikin and Schafer 
12] and, consistently, the RAMOD3 master equations, 
once converted into a coordinate form, recover the ana- 
lytical linearized case discussed by Kopeikin and others. 



Consequentely, those master equations are solvable with 
the same technique and may include any detectable rela- 
tivistic effects associated to the light propagation in that 
context, but only at e 2 level of accuracy. 



Instead, the RAMOD4 model, i.e. the case of a dy- 
namical space-time, fully preserves the "active" contents 
of gravity. Let us remind that formulas (|3T[) and (|38|) 
are not contemplated in other approaches: only equation 
(|3"5)) . referred to the spatial components, can be reduced, 
once transformed into a coordinate form, to the equation 
of propagation of photons in a weak gravitational field in 
the first pM approximation, as demonstrated in this pa- 
per. Actually, RAMOD4 master equations define a set of 
nonlinear differential equation for the four-vector £, the 
spatial local-line-of-sight as measured in the rest-space of 
a local barycentric observer which feels the gravitational 
field according to (ITU1) . This family of observers can de- 
fined in each event of the space-time, then also at the 
observer's location. Methods to solve light tracing an- 
alytically in RAMOD4 are under investigations; to this 
class of methods belongs the semi-anlytical solution of dc 
Felice and Preti I35B . 



A new analytical solution for the RAMOD master 
equations will surely represent a substantial improve- 
ment for the light propagation problem, since RAMOD 
"recipe", based on General Relativity, can tackle light 
tracing and detection to any order of accuracy and, even 
more important, consistently with the space-time geo- 
metrical set-up. This fact reflects the main purpose of the 
RAMOD approach which is to express the null geodesic, 
or the equivalent master equations, by all the physical 
quantities entering the process of observation, in order 
to entangle globally all the possible interactions of light 
with the background geometry. Keeping the physical as- 
pects of the problem guarantees the consistency of the 
measured physical effects with the intrinsic accuracy of 
the space-time. A preferred role is played in RAMOD 
by the local-line-of-sight, the main unknown of the mas- 
ter equations, which represents what locally the observer 
measures of the collected light. At the time of observation 
it represents the boundary condition to solve uniquely the 
light path. In a second step, once the light is traced back, 
any hidden subtle relativistic effects can be split by the 
global solution case by case, also according to the differ- 
ent level of accuracy needed, and finally quantified by a 
suitable coordinates system. 



As a concluding remark, by recovering results already 
presented in the literature, what this paper clearly estab- 
lishes is the full potential of the RAMOD construct: its 
general covariant formulation represents a well defined 
framework where any desired advancement in the light 
tracing problem and its subsequent detection as physical 
measurement can be contemplated. 
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Appendix A: Appendix 

Let us choose the coordinate system ^ adapted to 
the spatial hypersurfaces r = constant and designate 
two hypersurfaces T\ and r 2 . In order to compare proper 
length and proper time on them, we have to evaluate the 
line element, that is: 

ds 2 — g a pdx a dx lB — 

= g odx a dx° + 2g oi dx°dx t + g lj dx l dx j . (Al) 

We know that: 



so 



lido*,!*) = -5° a e^ 
u a (x°,x l ) = -g ^ 



= gipu p = g. l0 u° + g ik u k = 



(A2) 



9t0 



and 



g%ku 



-eV°- 



Substituting the latter expressions in (1A1[) . and adding 



and subtracting the same quantity 



j~of(dx ) , we get: 



ds 1 = (dx ) 2 



u l u 3 



0>j 



dx 1 



rdx' 



(da? 



rdx' 



(A3) 



Manipulating the condition u a u a = —1, we find: 



g a/ 3U u 



a„,/3 _ 



-goog 



gij^u 3 = — 1 



on 



giju l u 3 — 1 - gooff' 



u l 1 
Hence, denoting N % = — (shift factor) and N = ± — 



(lapse factor), equation (|A1|) becomes: 



ds 2 



-{Ndx ) 2 + g i:j [dx 1 - N l d^ [dx 1 - N j dx° 



(A4) 



In what follows we drop the tilde. Let k be the null 
tangent vector to the photon's geodesic. It can be ex- 
pressed as (|25lorf 



(Bl) 



k a = l a -(u\ k)u c 



Then, the geodesic equation can be derivate, step by step, 
as it follows: 



d[l a - (u | k)u a ] 
dX 



+ - ( u | k)vP][P - (u | k)ui] = o 



d\ y 1 ' dX \DX 1 
T a M [l^P - (u | k^v? - (u | h)Vu p + (u | fc) Vu 7 ] = 

= r -(«|A)^--« a [*^ T V r « j9 ] + 
dX dX 

T a M [fP - (u | k)(fu~< + Pu p ) + (u | fc)Vu 7 ] = 
where 

k^k T V T up = l p l T W T up - (u | k)l p iip, 

and the symbol Du/DX is the absolute derivative of 
u along the curve with parameter A and dots indicate 
derivative with respect to the parameter of the curve. In 
fact: 



and 



u^iip = 



/ j rv T u p = o 



because of condition u a u a = — 1. Let us trasform the 
affine parameter A of the geodesic into the parameter a 

da = —(u | k)dX. 

Therefore, the derivative of IB 1 1 results as: 

dl a _ d[-(u | k)l a ] da 



dX 



da dX ' 



i.e. 



dl a 

(u | fc) 2 -— + (u | k)T a (k T V T u (j k^) = 
da 

(u | k) 2 — - (u I k) 2 l a [(l T + u T )V T up¥] = 
da 

dl a 

(u | k) 2 - (u | k) 2 l a [(l T + u T )(l p + t/)V T M = 

da 

fffa 

(u I k) 2 -— - (u | k) 2 l a (l T PV T u + Pur). 
da 
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So the above geodesic equation becames: 



(u I k) 2 -— + (u I k) 2 ^- - (l<*+u a )[Pl T V T up + l%] 
da da 



(u | k) 2 T^[PP + {Fv? + Pu 13 ) + u'V] = 

Finally, dividing each member by (u \ k) 2 we find the 
geodesic equation (|3"2|) . 



Appendix C: The case of null expansion 

By imposing the free expansion condition (42j 

P?H V(a up) = Q pa i a /t = 0, (CI) 

where 

P«=8 a p +u p u a , 
P2 = 8Z+u a vP, 



!JaP = Va/3 + h a p, 



9oo 1 - h 



-(l + hoo), 



and 



d( a goo „ T 
Vaiip = e ¥ 'C( Q 3 /3 ) up) - T aB u r 



2ffo 



we find 

Q pa = P C ;P? V (a m 



{5 a p + u p u a )(6P + u a uP)(d {a gp )Q - r^ ff or) 



because of 



pa p p^9m um = 



p CT 29 



up) = 0. 



After replacing for T^g its expression with the metric 
coefficients 

^g Xr {d a g\p + d p g Xa - d x g a p) 
equation IC1I results term by term: 
P«pP(eVd ia gp) ) = ^ [d (p9a)0 + e 2 *g 0{p d mgry )o (C2) 
+ a 2v g 0(a d p) g m + e 4tp g 0{cr g p}0 dog o 



r p r a 1 a/3 "t 



p-a-ap^r- — jd g pa + e^d^g^o 



(C3) 



s 3<p go(pd<?)goo + ^ 5v g 0{ry g p)0 d g 00 . 



Finally, by summing terms (|C3I) and (|C3[) . we obtain: 

1 , 1 r 

<d P a = -^dogpa + e v go( P d\o\g<j)o + ^ e v go(ag P )odogoo- 

It is easy to see that for p = i and a = j the last equation 
becomes: 

®ij = ^dogij, 



implying for condition IC II 

dogtj = 0. 
Instead, for p = and a = i: 



(C4) 



Ooi = ^e v d g i + ^e 3v {g od gio + goid goo) (C5) 

+ ^e 5 ^ {googio + goigoo) d goo 

1 1 

= ^dogoi + 2 e v good goi 

= ~e v d a g 0i (l + e 2 *g O0 )=0 



=> do9oi = 
And, finally, for p = and a = 

e o = e^dogoo 



(C6) 



^+e 2 ^oo + ie^( 5o o) 2 



1 2 

= ^do^oo (1 + goo) 
by which, again, from equation IC II we get: 
=> <9o.9oo = 0. 



(C7) 



Appendix D: Calculations for formula (|35[) 

The spatial null geodesic with respect to the expansion- 
free congruence Ca is: 

— + — - (l a + u a )(J p up) + T'SJlP + u?)(F + xP) = 
da da ' ' 

To the first order in /i we have: 
1. 



"n 



5?fl + -/ir, 
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2. 



du a 
da 



F<5 Q (1 + ±hoo)(dih o)(l + h 00 ) « S^(l%h o); 



3. 



u/3 = u T V T up w (5q ( 1 + 5 ' l oo ) V T u /3 



r?.?f'+r? r(i+^oo)+r^F(i+i/ioo)+r? (i+ftoo) = o 

where we have, at the first order in h: 



and 



h Qi = 



dohij = 0. 



9/30 



25oo\/ — Soo 



d T (-gao) 



1 « t^A ffAO 



1. \ rl 



-(l + /l o)(?7 / 30+%o)(9r/lOo)+9r/l/30-r^ / 377AO 



Now, for a = 0: 



-Pdihoo 



--l j (djh o) - -d h 00 = 



hoo + d T hf3 Q - T^prixo + 0{h 2 ) 



d o h O0 = 0, 



and for a = k (k = 1, 2, 3): 



up = u T \7 T up « -r)pod hoa + d Q h po + T^. 



—-V 

da 



1 

+ 2 7 ' 



kX 



(dihx+djhxi-dxhiMPPy 



Adding all the above terms, the spatial null geodesic be- 
comes, to the first order in h: 



(^ft. ao + d Q h\i - d\h l0 )l l + (dohxj + djh xo - dxh Qj )l 3 + 



f + 1 -(l%h 00 )5 S 



^(■xVpodohoo + dohpo + T%p) 



(Dl) 



-01 



l 1 + 5 7 a (l + r>o) 



= 



i.e. 



dl a 1 - r- 1 i r- 

j^ + -{l % dih m )8%- [l a +5% (l+-h o)J [V(d h i0 +r 0i ) 



d h\ + d h\ - d\h 



0. 



The terms like dihko are null because in our system of 
coordinates g a o = and d^hki — for equation (IC4|) . 
So, finally: 



(D2) 



+ -V^idihxj + d,h Xl - dxhijWl 3 ) - ~r) KA d x h 00 = 



feA i 
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